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Abstract: An updated analysis of all available neutrino oscillation evidence in 
Solar experiments including the latest SNOES, CC and NC data (254d live time, 
NaCL enhanced efficiency) is presented. We obtain, for the fraction of active os- 



cillating neutrinos: sin 2 a = (® NC - $cc)/(®SSM ~ §cc) = O^Ol^o nearly 20cx 



from the pure sterile oscillation case. The fraction of oscillating sterile neutrinos 
cos 2 a < 0.12 (ler CL). At face value, these results might slightly favour the exis- 
tence of a small sterile oscillating sector. In the framework of two active neutrino 
oscillations we determine individual neutrino mixing parameters and their errors and 
we obtain Am 2 = 7.01 ± 0.08 x 10~ 5 eV 2 , tan 2 # = 0.42+°;J 2 . The main difference 
with previous analyses is a better resolution in parameter space. In particular the 
secondary region at larger mass differences (LMAII) is now excluded at 95% CL. 
The combined analysis of solar and Kamland data concludes that maximal mixing is 
not favoured at ~ 4 — 5a. This is not supported by the antineutrino reactor results 
alone. We also estimate the individual elements of the two neutrino mass matrix, 
writing M 2 = m 2 I + Mq, we obtain (la errors): 

M -1U eV I +Q.29 4fi n+0.56 

V' 3 - io -0.35 ^- OU -0.44 
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1. Introduction 

The Sudbury Neutrino Observatory (SNO) collaboration has recently presented data 
relative to the NaCl phase of the experiment [1]. The addition of NaCl to a pure 
D 2 detection medium increases the detector's sensitivity to the neutral-current 
(NC) reactions within its fiducial volume, enhancing the NC detection efficiency of 
about a factor three with respect to the previous 'no-salt' phase. This new result 
adds to a list of successive, compatible, ever more accurate results. The year 2002 
was fundamental for neutrino physics and many long standing puzzles found their 
final solution. The data published in springtime by SuperKamiokande [2, 3] and SNO 
[4, 5] and the first results from the reactor experiment KamLAND [6] in December, 
confirmed the evidence accumulated in about 30 years of solar neutrino experiments, 
proved in a crystal clear way that neutrinos are massive and oscillating particles 
and selected the so-called Large Mixing Angle (LMA) between the different possible 
solutions of the Solar Neutrino Problem (SNP). 

After the publication of the first KamLAND data a new era has started, in 
which the main issue is no longer to prove the validity of the oscillation hypothesis, 
but rather to determine with the best possible accuracy the oscillation parameters 
entering the PMNS mixing matrix. As it has been pointed out in previous papers 
(see for instance [7-10]), the accurate knowledge of the "solar mixing parameters" 
(6*i2 and Am^ in the case of normal hierarchy) is essential not only for a better 
comprehension of this sector of neutrino physics, but also for the future study of the 
other mixing parameters. 
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In fact the neutrino parameters drive the subdominant effects in the long baseline 
and neutrino factory experiments and a lack of information about their exact values 
could seriously reduce the accuracy in the determination of the other elements of the 
mixing matrix and eventually of the CP violation parameter 5. Hence, improving the 
present accuracy in the determination of the solar mixing parameters is mandatory 
not only to understand the physical potential of future experiments but also in order 
to optimize their set-up in the best possible way. 

From this point of view, the situation was still quite unsatisfactory after the pub- 
lication of the first KamLAND data. The many different analysis of these data [11- 
13] agree on the fact that these results, together with all the solar neutrino exper- 
iments, essentially select two distinct regions in the solar mixing parameter plan, 
both corresponding to the LMA solution. These two regions are usually denoted as 
"low LMA" and "high LMA" . The mixing angle was not maximal, even if it was 
large and it was constrained in the interval 0.27 < tan 2 6 , 12 < 0.88 at the 3a level. 
One can see that the uncertainty on the value of the mixing angle was still quite 
significant. Moreover it was not possible to discriminate between the two different 
possibilities for the mass difference. It is clear that a step further in the direction of 
pinning down more accurately the mixing parameters was needed. 

The SNO experiment measures the 8 B Solar neutrinos via the reactions [15-18]: 
1) Charged Current (CC): u e + d — > 2p + e~, 2) Elastic Scattering (ES): v x + e _ — > 
v x + e~. 3) Neutral Current (NC): u x + d — > p + n + is x . The first reaction is sensitive 
exclusively to electron neutrinos. The second, the same used in SuperKamiokande 
(SK), is instead sensitive, with different efficiencies, to all flavors. Finally the NC 
reaction is equally sensitive to all active neutrino species. Hence SNO can measure 
simultaneously the electron and active non-electron neutrino component of the solar 
flux at high energies ( > 5 MeV). The novelty of the 'salt phase' is that it is now 
possible to distinguish clearly the NC events from the ES and CC ones and therefore 
to analyse the data without making use of the no-spectrum-distortion hypothesis 
[1,19]. 

The non-electron component is found to be ~ 5a greater than zero, the standard 
prediction, thus providing the strongest evidence so far for flavour oscillation in the 
neutral lepton sector: the agreement of the total flux, provided by the NC measure- 
ment with the expectations implies as a by-product the confirmation of the validity 
of the SSM [20-22]. 

In this work we present an up-to-date analysis of all available Solar neutrino evi- 
dence, including the latest global SNO results [1] (together with the data of previous 
SNO phases [4,5,23] and of the other experiments [2,24,25], and of the KamLAND 
data [6]. A similar analysis using in addition the NaCl enhanced SNO spectrum will 
be presented elsewhere [14]. As a result, the new SNO data make the discrimination 
between the two different Am 2 regions possible. In particular the secondary region 
at larger mass differences secondary region at larger mass differences (LMAII or high 
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LMA) is now excluded at 95% CL. The combined analysis of solar and KamLAND 
data continues to conclude that maximal mixing is not favored at an even larger 
significance level as before. Finally, we show that the achieved resolution in neutrino 
parameter space is good enough for an estimation of the individual elements of two 
neutrino mass matrix and their errors. As it will be shown below we obtain a square 
mass matrix: 

/ 9 na+0.29 n -, r+0.29 

m - iu ev I +Q.29 4fin +o.56 

This result can be sharpened and a mass matrix can already be given if a concrete 
value for the absolute neutrino mass scale is assumed. 

The structure of the present work is the following. First (section 2) we update 
some model independent results which put in a quantitative basis the extent of the 
deviations with respect to the standard non-oscillating case and the relative impor- 
tance of active/sterile oscillations. After Section 3 dedicated to general description 
of methods, we determine (Section 4) the allowed areas in parameter space in the 
framework of active two neutrino oscillations from a standard statistical analysis. 
Individual values for Am 2 and tan 2 9 with error estimation are obtained from the 
analysis of marginal likelihoods. 

2. Some Model Independent Results 

Different quantities can be defined in order to make the evidence for disappearance 
and appearance of the neutrino flavours explicit. From the three fluxes measured by 
SNO is possible to define two useful ratios, deviations of these ratios with respect 
to their standard value are powerful tests for occurrence of new physics. Here we 
update the computations of Ref.[26] for the values for §cc/®es an d $cc/®nc being 
extremely careful with the treatment of the correlations on the incertitudes. The 
inclusion or not of these correlations can affect significantly the results for these 
ratios (see table II in Ref. [1] for a complete list of systematical errors). The results 
we obtain from the new SNO data are similar to the old ones except for a strong 
decrease in the error bars in some cases. From the value from SNO rates [1] we obtain 



$cc 

— u.oyi_ 096 , 



fl «Q1 +0.150 

U.UJ±_q QQg, 



'ES 

a value which is ~ 2.1 a away from the no-oscillation expectation value of one. The 
ratio of CC and NC fluxes gives the fraction of electron neutrinos remaining in the 
solar neutrino beam at detection point. We obtain 

$cc = 0.305+ ' 030 



-0.0245 

^NC 

this value is nominally many standard deviations (~ 20<r) away from the standard 
model case [27]. 
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Finally, if in addition to SNO data we consider the flux predicted by the solar 
standard model one can obtain sin 2 a, the fraction of active oscillating neutrinos, 
again using the SNO data and fully applying systematic correlations, we find: 



rin»„ = »" c -* cc = 0.940313 (") 

^SSM — ^CC 

cos 2 a < 0.12 (la), (2.2) 



where the fraction of oscillating sterile neutrinos cos 2 a = 1 — sin 2 a. The SSM flux is 
taken as the 8 B flux predicted in the revised Ref. [21]. Although slightly increasing, 
the central value is still clearly below one (only-active oscillations case). Although 
electron neutrinos are still allowed to oscillate into sterile neutrinos the hypothesis of 
transitions to only sterile neutrinos is rejected at nearly 15a. On the other hand, as 
a consequence of the reduced error bars, this data can be taken as a mildly positive 
hint in favour of a small sterile component: the pure active case is now one sigma 
away from the central value (pure active oscillations are "excluded" at one sigma). 



3. Methods and statistical procedures 

The computation of the neutrino oscillation probabilities in Solar and Earth matter 
and of the expected signal in each experiment follows the standard methods found in 
the literature [28-36]. We solve numerically [29, 37], the neutrino evolution equations 
for all the oscillation parameter space. The survival probabilities for an electron neu- 
trino, produced in the Sun, to arrive at the Earth are calculated in three steps. The 
propagation from the production point to sun's surface is computed numerically in 
all the parameter range using the electron number density n e given by the BPB2001 
model [21] averaging over the production point. The propagation in vacuum from 
the Sun surface to the Earth is computed analytically. The averaging over the annual 
variation of the orbit is also exactly performed using simple Bessel functions. To take 
the Earth matter effects into account, we adopt a spherical model of the Earth density 
and chemical composition. In this model, the Earth is divided in eleven radial density 
zones [38], in each of which a polynomial interpolation is used to obtain the electron 
density. The composition of the neutrino propagation in the three different regions 
is performed exactly using an evolution operator formalism [36]. The final survival 
probabilities are obtained from the corresponding (non-pure) density matrices built 
from the evolution operators in each of these three regions. The night quantities are 
obtained using appropriate weights which depend on the neutrino impact parameter 
and the sagitta distance from neutrino trajectory to the Earth center, for each de- 
tector's geographical location. In this analysis in addition to night probabilities we 
will need the partial night probabilities corresponding to the 6 zenith angle bin data 
presented by SK [3]. They are obtained using appropriate weights which depend on 



-4- 



the neutrino impact parameter and the sagitta distance from neutrino trajectory to 
the Earth's center, for each detector's geographical location. 

The expected signal in each detector is obtained by convoluting neutrino fluxes, 
oscillation probabilities, neutrino cross sections and detector energy response func- 
tions. We have used neutrino-electron elastic cross sections which include radiative 
corrections [39]. Neutrino cross sections on deuterium needed for the computation 
of the SNO measurements are taken from [40]. Detector effects are summarized by 
the respective response functions, obtained by taking into account both the energy 
resolution and the detector efficiency. We obtained the energy resolution function 
for SK using the data presented in [42-44] . The effective threshold efficiencies, which 
take into account the live time for each experimental period, are incorporated into 
our simulation program. They are obtained from [45]. The resolution function and 
other characteristics for SNO used here are those given in Refs.fl, 4, 5, 41]. 

The two principal ingredients in the calculation of the expected signal in Kam- 
LAND are the reactor flux and the antineutrino cross section on protons. We refer 
to Ref. [11] for a detailed account of the methods used here. In summary, in order to 
obtain the expected number of events at KamLAND, we sum the expectations for all 
the relevant reactor sources weighting each source by its power and distance to the 
detector (table II in Ref. [46] ), assuming the same spectrum originated from each 
reactor. We sum over the nearby power reactors, we neglect farther Japanese and 
Korean reactors and even farther rest-of-the-world reactors which give only a minor 
additional contribution. The average expected signal in each energy bin is given by 
the convolution of the oscillation probability averaged over the distance and power 
of the different reactors. Expressions for the antineutrino capture cross section are 
taken from the literature [47,48]. The matrix element for this cross section can be 
written in terms of the neutron half-life, we have used the latest published value 
ti/2 — 613.9 ±0.55 [50]. The antineutrino flux spectrum, the relative reactor-reactor 
power normalization which is included in the definition of the effective probability 
and the energy resolution of KamLAND are used in addition. The energy resolution 
in the prompt positron detection is obtained by us from the raw calibration data 
presented in Ref. [11, 51]. Moreover, we assume a 408 ton fiducial mass and standard 
nuclear plant power and fuel schedule, we take an averaged, time- independent, fuel 
composition equal for each detector. Detection efficiency is taken close 100% and 
independent of the energy [6]. 

3.1 Statistical Analysis 

The statistical significance of the neutrino oscillation hypothesis is tested with a 
standard x 2 method which is explained in detail in Ref. [29]. Our present analysis 
is based on the consideration of a x 2 sum of two distinct contributions, one coming 
from SNO global rates and all the rest of solar neutrino data and the contribution 
of the KamLAND experiment x 2 = X© + Xkl- The contribution from KamLAND 
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includes the binned signal (See table 2 in Ref. [11]) as is explained in detail in Ref. [6]. 
In summary, the KamLAND contribution is made of two parts, one with x 2 g i kl = 
(R exp — R th ) 2 /a 2 . The experimental signal R exp and statistical and systematic errors 
are shown in Table [I]. The contribution of the KamLAND spectrum is as follows: 

xl pec ,KL = («# th - # CXP ) 4 Knc + o- 2 corr )~ l (aR th - R C ^) (3.1) 

The total error matrix a is computed as a sum of assumed systematic deviations, 
&sys/S ~ 6.5%, mainly coming from flux uncertainty (3%), energy calibration and 
threshold (see table II of Ref. [6]. for a total systematic error ~ 6.4%), see also 
Ref. [46, 51, 56]) and statistical errors. The parameter a is a free normalization pa- 
rameter. The effect of systematic sources on individual bin deviations has been 
computed by us studying the influence on the response function, furtherly we have 
assumed full correlation among bins. 

The solar neutrino contribution can be written in the following way: 

X© = Xgiob + X 2 sk + Xsno- ( 3 - 2 ) 

The function Xgiob correspond to the total event rates measured at the Homes- 
take experiment [24] and at the gallium experiments SAGE [52], GNO [53] and 
GALLEX [54]. We follow closely the definition used in previous works (see Ref. [26] 
for definitions and Table (1) in Ref. [26] for an explicit list of results and other ref- 
erences). The contribution to the x 2 from the SuperKamiokande data (Xsk) nas 
been obtained by using double-binned data in energy and zenith angle (see table 2 
in Ref. [3] and also Ref. [2]): 8 energy bins of variable width and 7 zenith angle bins 
which include the day bin and 6 night ones (see Ref. [6]). 

The contribution of SNO to the x 2 is given by xIno = X 2 g i,sNO + X 2 sp ec,SNO where 
X%ec,SNO i s the spectrum contribution made up by the day and night values of the 
total (NC+CC+ES) SNO signal for the different values of the spectrum [11]. The 
new component corresponding to the individual global signals is given by 

4.SNO = E ( aRth - ^ Kat + ^syst)" 1 (*R th ~ i* CXP ), (3.3) 

i=ES,CC,NC 

where the signal R vectors of dimension 3 are made up by the values of the total 
ES, CC, NC SNO signals. The statistical contribution to the covariance matrix, 
cr sta t is non-diagonal since the different fluxes are derived from a fit to a single data 
sample [1,41]. The part of the matrix related to the systematical errors includes 
contributions from neutron capture efficiency and other geometrical inefficiencies 
appearing in the statistical separation of ES, CC and NC events as presented in 
Ref. [41]. 
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4. Results and Discussion 

To test a particular oscillation hypothesis against the parameters of the best fit 
and obtain allowed regions in parameter space we perform a minimisation of the 
three dimensional function x 2 (Am 2 , tan 2 9, a). For a = a m i n , a given point in the 
oscillation parameter space is allowed if the globally subtracted quantity fulfills the 
condition Ax 2 = x 2 (AmV) - xLin < xl{CL). Where xjU 3 (90%, 95%, ...) are the 
quantiles for three degrees of freedom. 

The results are shown in Figs.|T] where we have generated acceptance contours in 
the Am 2 -tan 2 9 plane. In Fig.^|we present the same results as in Fig.|l]but using linear 
scales. The resolution in the neutrino parameter space has become good enough for 
this to become useful. With the actual experimental precision we can assert that 
the mixing parameters are now known much better than as order of magnitude only. 
In Table (0) we present the best fit parameters or local minimum obtained from the 
minimisation of the full \ 2 function. 

The main difference with previous analysis is a better resolution in parameter 
space. The previously two well separated solutions LMAI,LMAII have now disap- 
peared. In particular the secondary region at larger mass differences (LMAII) is now 
excluded at 95% CL. 

The introduction of the new solar data in general strongly diminishes the favored 
value for the mixing angle with respect to the KamLAND result alone [11]. The 
final value is more near to those values favored by the solar data alone than to 
the KamLAND ones. As an important consequence, the combined analysis of solar 
and KamLAND data concludes that maximal mixing is not favored at ~ 4 — 5a. 
This conclusion is not supported by the antineutrino, earth-controlled, conceptually 
simpler KamLAND results alone. As we already pointed out in Ref.[ll], this effect 
could be simply due to the present low KamLAND statistics or, more worrying, 
to some statistical artifact derived from the complexity of the analysis and of the 
heterogeneity of binned data involved. 

Additionally, we perform a second kind of analysis in order to obtain concrete 
values for the individual oscillation parameters and estimates for their uncertainties. 
We study the marginalised parameter constraints where the x 2 quantity is converted 
into likelihood using the expression C/Co = e~( x2 ~ x ™ in ^ 2 . This normalized marginal 
likelihood, obtained from the integration of C for each of the variables, is plotted 
in Figs. ([|) for each of the oscillation parameters Am 2 and tan 2 9. Concrete values 
for the parameters are extracted by fitting one- or two-sided Gaussian distributions 
to any of the peaks (fits not showed in the plots). In both cases, for angle and the 
mass difference distributions the goodness of fit of the Gaussian fit to each individual 
peak is excellent (g.o.f ~ 100%). The values for the parameters obtained in this way 
appear in Table 0. The errors obtained from this method are assigned to the \ 2 
minimisation values. The central values are fully consistent and very similar to the 
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values obtained from simple \ 2 minimisation. In particular, the maximal mixing 
solution is again excluded at the ~ 4 — 5cr level. A common feature to previous 
analysis presented by us [11] is that, although both are mutually compatible, the 
slight difference of the value obtained for the mixing angle is well explained by the 
shape of the allowed regions in Fig [I] (right): the right elongation of these shift the 
value of the integral which defines the marginal distribution for tan 2 9. Additional 
variability can be easily introduced if would have used different prior information or 
mixing parameterizations. 

We will again use the technique of marginal distributions in the next paragraphs 
to obtain an estimation of the individual elements of the neutrino mass matrix and 
their errors. 

4.1 An estimation of the neutrino mass matrix 

The square of the neutrino mass matrix can be written in the flavour basis as M 2 = 
UMpU^ where Md is diagonal and U is an unitary (purely active oscillations are 
assumed) mixing matrix. Subtracting one of the diagonal entries we have 

M 2 = mil + M 2 = mil + UM%U\ 

where I is the identity matrix. In this way we distinguish in the mass matrix a part, 
Mg, which affects and can be determined by oscillation experiments and another 
one, mil, which does not. Evidently, the off-diagonal elements of the mass matrix 
are fully measurable by oscillation experiments. 

Restricting ourselves for the sake of simplicity to two neutrino oscillations, we 
have 

,^9 9 T ,,9 9 T 9 / sin 2 9 sin 9 cos 9 \ , „ „ . 

M 2 = mil + M 2 =mll + Am 2 ( 4.1 

\sm0cos0 cos z 9 / 

with Am 2 = m\ — m\. The individual elements of the matrix M can simply be 
estimated from the oscillation parameters obtained before. For example for tan 2 9 ~ 
0.40, Am 2 ~7x 10" 5 eV 2 we would obtain (M 2 ) 22 ~5x 10" 5 eV 2 . 

Our objective is however to estimate how well the individual errors of the mass 
matrix can be extracted already at present by the existing experimental evidence. 
For this purpose we have applied similar arguments as those used before to obtain 
marginal distributions and errors for individual parameters from them. Using again 
as likelihood function the quantity C/Co(Am 2 , tan 2 9) = e~ ( - x ' 2 ~ x ™ i ™^ 2 we obtained 
the individual probability distributions for any of the elements of the matrix M Q . 
Average values and la errors are obtained from two-sided Gaussian fits to these 
distributions. 

From this procedure we obtain: 

/ 2f)fi +0.29 O1P-+0.29X 



- 8 - 



One can go further supposing a concrete value for m 2 from elsewhere. If we take 
m\ » Am 2 then we can directly write the mass matrix 

M = mj + ^—Ml (4.3) 



Supposing for example m\ = \ eV 



M ~ V i.6ol8:lf io- 5 i.oo + 4.601°^ io- 5 J ■ l4 - 4j 

this is the final two neutrino mass matrix which can be obtained from present oscil- 
lation evidence coming from solar and reactor neutrinos. 



5. Summary and Conclusions 

In this work we have presented an up-to-date analysis of all available Solar neutrino 
evidence including latest SNO results with NaCl enhanced efficiency in the most sim- 
ple framework. The increasingly accurate direct measurement via the NC reaction 
on deuterium of 8 B neutrinos combined with the CC results have largely confirmed 
the neutrino oscillation hypothesis. 

In a model independent basis, We obtain the following values for the ratios: 
§cc/$nc = 0.30518:^2, ®cc/®es = 0.6911^96- The fraction a of oscillating neu- 
trinos into active and sterile ones are computed to be: 

sin 2 a = 0.940ig;|jaj, cos 2 a < 0.12 (la), (5.1) 

where the fraction of oscillating sterile neutrinos cos 2 a = 1— sin 2 a. Although slightly 
increasing, the central value is still clearly below one, the only-active oscillations 
case. The hypothesis of transitions to only sterile neutrinos is well rejected but, as 
a consequence of the reduced error bars, this data can be taken as a mildly positive 
hint in favour of a small sterile component: the pure active case is now one sigma 
away from the central value. 

We have obtained the allowed area in parameter space and individual values for 
Am 2 and tan 2 9 with error estimation from the analysis of marginal likelihoods. We 
have shown that it is already possible to determine at present active two neutrino 
oscillation parameters with relatively good accuracy. In the framework of two active 
neutrino oscillations we obtain 

Am 2 = 7.01 ± 0.08 x 10" 5 eV 2 , tan 2 9 = 0.42l°;J 2 . 

The combined analysis of solar and KamLAND data concludes that maximal mixing 
is not favored at ~ 4 — 5cr. This conclusion is not supported by the antineutrino, 
earth-controlled, conceptually simpler KamLAND results alone. 
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We estimate the individual elements of the two neutrino mass matrix, we show 
that individual elements of this matrix can be determined with an error ~ 10% from 
present experimental evidence. 
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New SNO data[lj: 






SNO-ES 


2.21 ±0.28 ±0.10 


0.406 ±0.091 


SNO-CC 


1.59 ±0.08 ±0.07 


0.292 ± 0.056 


SNO-NC 


5.21 ± 0.27 ± 0.38 


0.958 ±0.193 


Other Solar data: 






D 2 SNO-ES [4, 5] 


2.39 ±0.23 ±0.12 


0.439 ±0.092 


D 2 SNO-CC [4, 5] 


1.76 ±0.05 ±0.09 


0.324 ±0.054 


D 2 SNO-NC [4, 5] 


5.09 ±0.43 ±0.45 


0.936 ±0.208 


SK [2] 


2.32 ±0.03 ±0.08 


0.451 ±0.011 


CI [55] 


2.56 ±0.16 ±0.16 


0.332 ±0.056 


SAGE [52] 


67.2 ±7.0 ±3.2 


0.521 ±0.067 


GNO-GALLEX [53, 54] 


74.1 ±6.7 ±3.5 


0.600 ±0.067 



Table 1: Summary of data used in this work. The observed signal (Spata) and ratios 
SData/ Sssm with respect to the BPB2001 model are reported. The SK and SNO rates 
are in 10 6 cm~ 2 s _1 units. The CI , SAGE and GNO-GALLEX measurements are in SNU 
units. In this work we use the combined results of SAGE and GNO-GALLEX: Sq^/S^m 
(Ga = SAGE±GALLEX±GNO)=0.579 ± 0.050. The SSM 8 B total flux is taken from the 
(revised) BPB2001 model [21]: <M 8 B ) = 5.44(1+°;^) x 10 6 cm" 2 s _1 . In addition we 
have used from reactor Kamland measurements the signal ratio R = 0.611 ± 0.085 ± 0.041 
[6] and its signal spectrum [6, 11]. 



Am 2 (eV 2 ) tan 2 # 
From minimization x 2 '- 7.01 x 10 -5 0.42 
From Marg. Fit, ( ±la): 7.30±g;g| x 10" 5 0.46j[ 



Table 2: Mixing parameters: from \ 2 minimization, \ 2 /ndf = 0.94 (Fig. || right) and 
from double-sided fit to the peak of marginal likelihood distributions (Figs.||). 




Figure 1: Allowed areas in the two neutrino parameter space. The colored areas are the 
allowed regions at 90, 95, 99 and 99.7% CL relative to the absolute minimum. (Left) Solar 
evidence (CL,GA,SK,SNO,SNO-salt). (Right) Kamland spectrum plus solar evidence. 
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Figure 2: Kamland spectrum plus solar evidence as Fig.[T] (right), in linear scale which 
allows for a better comparison between the different regions. 
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Figure 3: Marginalized likelihood distributions for each of the oscillation parameters Am 2 
(right), tan 2 9 (left) corresponding to the solar plus Kamland evidence (Fig.[j](Right)). The 
curves are in arbitrary units with normalization to the maximum height. Values for the 
peak position are obtained by fitting two-sided Gaussian distributions (not showed in the 
plot). See Table 2 for values of the position and widths of the peaks. 
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